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Abstract 

The standard determination of the QED coupling on the Z pole is performed using 
the latest available data for R. The direct application of analytic continuation techniques 
5_i ■ is found not to improve the accuracy of the value of a(M|). However they help to 

resolve an ambiguity in the values of R in the energy region y/s < 2 GeV, which, in turn, 
reduces the uncertainty in a(M§). Moreover, they provide a sensitive determination of 
the mass of the charm quark. The favoured solution, which uses the inclusive data for R 
for yfe < 2 GeV, has a pole mass m c = 1.33 - 1.40 GeV and a _1 (M|) = 128.972 ± 0.026; 
whereas if the sum of the exclusive channels is used to determine R in this region, we find 
or l (Ml) = 128.941 ±0.029. 



1 Introduction 



The value of the QED coupling at the Z boson mass, a(M|), is the poorest known of the three 
parameters necessary to define the standard electroweak model, which, for example, may be 
taken to be Gf, Mz and a(M z ). The value of a(M§) is obtained from 

a' 1 = a(O)" 1 

using the relation 

= (l - Aai e p(s) 

where the leptonic contribution to the running of the a is known to 3 loops |l] 

Aa lcp (Mf) = 314.98 x 1(T 4 . (3) 

^From now on we omit the superscript (5) on A«had and assume that it corresponds to five 
flavours. We will include the contribution of the sixth flavour, Aa top (Mf ) = —0.76 x 10~ 4 , at 
the end. To determine the hadronic contribution it is traditional to evaluate 

as „ f°° R(s')ds' 

at s — Mz, where R = <r(e + e~ — > hadrons)/cr(e + e~ — > 

The main uncertainty in the calculation of Aahad comes from the lack of precise knowl- 
edge of R(s') in the energy region 1.5 < vV^S 3 GeV, see Fig. 1. In the upper half of this 
interval the situation has recently improved with the new (preliminary) BES-II measurements 
0. Nevertheless there remains a major problem due to the discrepancy between the inclusive 
measurements of e + e~ — > hadrons and the value of the cross section deduced from the sum of 
all the exclusive hadronic channels (e + e~ — > 2tc, 3tc, . . . , KK, . . .), see Fig. 1. 

Recently dispersion relation (f|) has been re-evaluated at s = M| || |J, incorporating the 
new BES-II data for R(s'). In Section 2 we give the details of the determination of Ref. and, 
in particular, expose the dilemma with the input values of R(s') in the region ^fs 1 < 2 GeV. 
In Section 3, following Jegerlehner ||, we describe an attempt to better determine Aa had (M|) 
by evaluating dispersion relation (f|) in the space-like region, at s = — s say, and then using 
perturbative QCD to analytically continue from s = —So —>■ — M§ — > M§. Although this 
procedure is found to reduce the error associated with the data for R(s'), it is more than 
compensated by the uncertainties in the analytic continuation coming from the choice of the 
mass of the charm quark and the QCD scale. 

However analytic continuation offers the possibility to resolve the dilemma in the data for 
R(s') in the region \fs 1 < 2 GeV (see Section 4), and to give a reasonably accurate determination 
of the pole mass m c of the charm quark (see Section 5). Clearly a resolution of the dilemma 
will improve the direct determination of Aahad(^l) obtained by evaluating ([|) at s = Mf . In 
Section 6 we present our conclusions. 
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Figure 1: The quantity R(s) versus yfs in the critical low energy interval, ^J~s < 3 GeV. The band 
below y/s = 2.125 GeV now illustrates the bounds of the summed exclusive channels. The inclusive 
data are explicitly plotted, and above ^fs = 1.46 GeV the curve shows the central value of their 
interpolation. In the overlapping interval there is a distinct discrepancy between the two (in principle) 
complementary measurements. The central perturbative QCD prediction at O(ag) is plotted through 
the inclusive region for comparison. Finally, the vertical lines denote the central positions of the </> 
and J/tp resonances. (See the note added in proof for the final BES measurements |pq].) 
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2 Direct determination of Aa^^M^) 



In this section we give the details of the recent determination^ of Aahad(^l) that was presented 
in P|- We evaluated dispersion relation (|) at s = Mf using the experimental data |2|, || [7|, ||. [i| 
for R(s') in the intervals 2m^ < \fs' < 2.8 GeV and 3.74 < yfs* < 5 GeV, together with the 
J/tp,ip' and T resonance contributions. In the remaining regions (2.8 < yfs' < 3.74 and 
\fs' > 5 GeV) we calculate R(s') from perturbative QCD using the two-loop expression with 
the m c and nib quark masses included and the massless three-loop expression [|1(J calculated 
in the MS renormalization scheme?]. We estimate the 'perturbative' error on R(s') by allowing 
m c , m b , M z to vary within the uncertainties quoted in JT^], by taking a,s(Mf ) = 0.119 ± 0.002 
and by varying the scale of cts{cs) in the range 0.25 < c < 4. 

The errors on the 'data' values of R{s') are calculated using a correlated x 2 minimization to 



combine the different data sets, as described in detail in Ref. ||13|| . The data, together with the 
error band used in the 3.74 < < 5 GeV interval, are shown in Fig. 2. For \/~P < 1.46 GeV 
the sum of the data for the exclusive channels is used to compute R(s'), see Table 1. Recently 
there have been improvements in our knowledge of the exclusive channels. This can be seen, for 
example, in the data |7J for the 2tt channel shown in Fig. 3, or the data || for the 4n channel 
shown in Fig. 4. 

For -v/s 7 > 1.46 GeV we also have inclusive measurements of R(s'). These differ significantly 
from the sum of the exclusive channels, see Fig. 1. This poses a dilemma. The new (preliminary) 
BES-II data 0, which extend down to yfP = 2 GeV, appear to match better to the inclusive 
measurements, but the distinction is not conclusive. We therefore, throughout this paper, 
take two alternative choices of the data in the interval 1.46 < \fs' < 1.9 GeV. We first use 
the inclusive data and then we repeat the analysis using the exclusive data (with the error 
band shown in Fig. 1). For simplicity, we refer to these as the 'inclusive' and 'exclusive' data 
choices. In the later sections of this paper we study ways to resolve this dilemma and we 
present evidence which favours the 'inclusive' behaviour of R(s') in this interval. In Table 2 
we list the contributions to the dispersion relation ([|) from specific y/s' intervals for both the 
above choices of data. In the Table we also include the corresponding values of Aa[^ d (M|) and 
a _1 (M|). We see that the ambiguity in the input for R(s') in the region y/s' < 2 GeV itself 
leads to an uncertainty of the size of the quoted errors on Aa^ d (M|). We attempt to resolve 
this ambiguity in Section 4. 

The values that we obtain for a _1 (M|) are compared with other recent determinations 
in Fig. 5. We also include on this plot two 1994- 5 determinations in order to gain some 
insight. First, we show the value obtained by Martin and Zeppenfeld [|15| which made use 
of perturbative QCD, as has become common practice, and which used 'inclusive' data for 
s' > 1.46 GeV and rescaled data in the charm resonance region. Second, we show the value 



1 A correction to the analysis of Ref. Q shifts the value of Aahad by 0.44 x 10 4 . 

2 The uncertainty due to using a different scheme may be estimated to be of the order of the 0(ag) correction, 
iT^eir^ias j V) 4 . We may take r% = —128 jllj which leads to a negligible uncertainty in R(s'). 
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Final state 


A «Lad( M i) • !o 4 

2m„ - 1.46 GeV 


1.46-1.9 GeV 




33.93 ± 0.52 


0.17± 0.06 




0.30 ± 0.04 


0.17 ± 0.05 


7T^~ 7T 7r^7T^ 


2 00 ± 08 


2 99 ± 31 


OJ 7T° « 


0.12 ± 0.02 


0.04 ± 0.01 




1.45 ±0.05 


2.29 ± 0.09 




0.09 ±0.04 


0.70 ±0.25 




0.04 ±0.05 


0.33 ±0.22 


0; 7T + 7T~ ^ 




0.02 ±0.00 


7T + 7r _ 7r + 7r _ 7r + 7T _ 




0.05 ±0.02 


7T + 7T - 7T + 7T ~~ 7T° 7T° 


0.02 ±0.01 


0.82 ±0.09 




0.01 ±0.01 


0.61 ±0.61 




0.02 ±0.02 


0.12 ±0.04 




0.53 ±0.05 


0.16 ±0.02 




0.15 ±0.11 


0.04 ±0.02 




0.03 ±0.01 


0.28 ±0.05 


K+K-tt 




0.10 ±0.07 


K° s K° L Tr° ( 3 ) 




0.10 ±0.07 


JK Ktttt W 


0.01 ±0.25 


1.04 ±0.67 


Sum of contributions 


38.76 ±0.79 


10.32 ± 1.06 



Table 1: A detailed breakdown of the individual exclusive channel contributions to Aa[^ d (M§). The 
dominant contribution arises from the e + e~ — > ir + ir~, and the next most significant contributions 
are obtained from e + e~ — > Tr + ir~iT + ir~ and e + e~ — ► 7r + 7r~7r°7r°, depicted in Fig. ^. The channels 
marked with W have been corrected for missing modes. The channel highlighted by ^ has had the 
77 — > 3tt contribution subtracted. Those modes marked by ( 3 ) have their contributions deduced from 
isospin relations. The modes described in ( 4 ) are deduced from the 'partially' inclusive measurements 
of e + e~ — > Kg + X, with modes explicitly included elsewhere subtracted. We have checked the 
contributions to the cross-section from each annihilation channel with the detailed decomposition 



given in [13|, and find excellent agreement between the two evaluations. 
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Figure 2: The quantity R(s) in the vicinity of the charm threshold 3.74 < i/s < 5 GeV. The Mark 
1, DASP and PLUTO data have been scaled by factors of 0.84, 0.88 and 0.95 so as to agree with 
the perturbative QCD prediction in the continuum regions safely above and beneath threshold. To 
guide the eye, vertical lines denoting the position of the -0(4040), -0(4160) and -0(4415) resonance 
centres have been superimposed. The band illustrates the interpolation derived from the compilation 
of the (rescaled) data. The perturbative prediction for R to 0(a^) is depicted in the continuum. The 
evaluations of Aahad in this work use the perturbative prediction for R(s) in the regions 2.8 < yfs < 
3.74 GeV and ^/s > 5 GeV. 
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Figure 3: The cross-section for pion pair production, a n7T (s), versus around the p-resonance region, 
2m n < \fs < 1 GeV. The data @ include the recent, accurate results from Novosibirsk. The band 
illustrates the spread of uncertainty about a central value interpolated from the data compilation. 



The line at low energies shows the chiral expansion of the two pion cross-section [14]. 
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of Eidelman and Jegerlehner Jl6[ which was obtained using data in all intervals, and hence the 
larger errors. For interest, we compare the individual contributions and errors of our present 
'inclusive' determination with those of the 1995 analysis of Eidelman and Jegerlehner in Table 3. 

In Fig. 6 we show the y 2 profiles^ obtained using the 'inclusive' and 'exclusive' determi- 
nations of the QED coupling a(M§) in fits to the latest compilation of electroweak data for 
different values of the mass of the (standard Model) Higgs boson. We see that the minimum 
obtained using the 'inclusive' value, a(M§) = 1/128.972, is close to the LEP2 bound on the 
Higgs mass. 



3 Analytic continuation in the space-like region 



There have been several studies [17], [j| of analytic behaviour in the complex s-plane in attempts 
to reduce the dependence of the determination of Aahad(^l) on the observed values of R in 
the region in which it is poorly known. These techniques have been reviewed by Jegerlehner ||. 
He concludes that it is difficult to reduce the error on Aa had due to the data in this way. He 
advocates the following analytic continuation method to determine Aahad(^I)- First, evaluate 
(fl) for space-like s = —sq and then use perturbative QCD to continue to s = — M|, that is 



Aa h . 



ad I 



-M 2 Z ) 



Aa had (-M|) - Aa had (-s ; 



QCD 



+ Act had (-s ) 



data. 



(5) 



3 We thank Martin Griinewald for making this plot. 
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y/s interval ( GeV) 


Aa had( M I) ' 104 contribution 


Origin of contribution 


2m n - 1.46° 


»■■»*{%} 


Pion form factor data 


1.46 - 1.90 


\ 


f 8.62±0.60 c 
10.32 ± 1.06 


J Inclusive data 
\ Exclusive summation 


1.90 - 2.80 


\ 


' 13.26 ±0.83 c 
13.79 ±0.83 


J Inclusive data 
Exclusive summation 


2 80 - 3 74 




9.73±0.05 d 


Perturbative QCD 


3.74 - 5.00 




15.13 ±0.36 


Charm data 


5.00 - oo 




169.97 ±0.64 d 


Perturbative QCD 


uj, <j>, tp's, Y's 




18.79 ±0.58 


Brcit-Wigner resonances 


a4 5 2j(m1).io 4 


1 


' 274.26 ± 1.90 
276.49 ±2.14 


j Inclusive data 
\ Exclusive summation 


a-\Ml) 


f 128.972 ±0.026 
\ 128.941 ±0.029 


( Inclusive data 
[ Exclusive summation 



Table 2: The individual contributions to the hadronic component of the shift in fine structure constant, 
Aa{^ d (M|) • 10 4 . The upper (lower) error in the result labelled a corresponds to the 2ir (remaining) 
exclusive channels. Contributions labelled with superscripts b , c and d have common error sources 
which are added linearly. Remaining errors are added in quadrature. 



where So is chosen sufficiently large (y/so > 2 GeV) for the QCD contribution in square brackets 
to be known accurately^, such that the error in Aah a d( — dominantly reflects the error in 
the data for R(s'). The error associated with the final continuation round the semicircle to 
A« had (M|) is negligible 

Aa had (Mf) = Aa had (-M|) + (0.42 ± 0.02) x 10" 4 . (6) 

Jegerlehner chose y/so = 2.5 GeV and foundf] 

Aa had (M|) = (277.82 ± 2.54) x 10" 4 (7) 

where the error was entirely attributed to that for the contribution Aahad( — so) data to (|5p. 

We will examine this proposal below. In particular we will investigate whether it is possible 
to develop this technique either to select between the inclusive/exclusive R(s') data choices 
in the region < 2 GeV, or to reduce the importance of the data contribution (and its 
associated error) from this domain. 

Suppose, for example, we evaluate a(Mf ) from ([|), @ and (H for a range of different values 
of Sq. In principle, we should always get the same answer. If the answer varies significantly 
either the data for R(s') is not quite correct or the theory input is deficient in some way or, 

4 Previous studies 01| had indicated how large sq had to be to avoid uncertainties due to parton condensate 
contributions. 

5 The recent BES-II data [|| were not available for the analysis of Ref. ||. 
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Final state 


y>s interval (GeV) 


Contribution trom |16| 


Current evaluation 


P 


2/71^ - 0.81 


26.08 ± 0.68 


25.32 ± 0.52 


w 


0.42 - 0.81 


2.93 ± 0.09 


3.07 ± 0.10 




1.00 - 1.04 


5.08 ± 0.14 


5.08 ± 0.19 


JH> 




11.34 ± 0.82 


9.41 ± 0.53 (+1.93=11.35) 


T 




1.18 ± 0.08 


1.22 ± 0.04 


hadrons 


0.81 - 1.40 


13.83 ± 0.80 


12.24 ± 0.54 








r 1.40- 1.46 1.21 ±0.07 Exc. \ 


hadrons 


1.40 - 3.10 


27.62 ± 4.02 


\ 1.46- 2.8 21.88 ±1.43 Inc. > 








1 2.8 - 3.10 3.43 ±0.02 pQCD. J 


hadrons 


O 1 A o /in 

3.10 - 3.60 


5.82 ± 1.16 


5.02 ± 0.03 








( % Rn q 7/1 i oa x n m r,npn "\ 

1 O.OU — 0./4 l.zo ± U.U1 p^^U. I 


hadrons 


3.60 - 9.46 


50.60 ± 3.33 


< 3.74-5.0 15.13 ±0.36 Inc. 1 








1 5.0 - 9.46 35.51 ±0.21 pQCD. J 


hadrons 


9.46 - 40.0 


93.07 ± 3.50 


91.77 ± 0.19 


perturbative QCD 


40.0 - oo 


42.82 ± 0.10 


42.70 ± 0.24 


Total 


2771^ - OO 


280.37 ± 6.54 


274.26 ± 1.90 



Table 3: A comparison of the individual contributions to Aa{^ d (M|) • 10 4 found in the 1995 'data- 



driven' analysis of Eidelman and Jegerlehner 
according to the energy intervals used in [16]. 



with those of our inclusive analysis, decomposed 



more likely, it is a combination of both. The interplay between the uncertainties in the theory 
and the data (that is, in the two terms on the right-hand-side of (|5])) play a crucial role in 
this type of analysis. If it is possible to find a choice of input data, together with a physically 
meaningful set of theory parameters (charm mass m c , choice of scale etc.), which give a stable 
value of a(M§) for different choices of so, then it will be a powerful argument in favour of their 
veracity. 

Indeed, imagine one extreme in which the theory contribution to ([5]) was known precisely; 
that is, there is no error associated with the term in square brackets. Then the behaviour of 
the variation of a(M|) as a function of Sq would highlight the domain (or domains) in which 
the data were wrong and, moreover, specify the approximate corrections that are necessary. 

In this section we evaluate Aa^adis) of (^) in the space-like domain s = — s (with s > 0) 
for a range of different values of Sq. For each s we then use perturbative QCD to perform 
the analytic continuation to s = M|, as given in (f|) and @. A sample of the results for 
A«had(— s o) is presented in Table 4, together with the conventional time-like evaluation of @ 
at s = +Mf . We see that the error on the space-like evaluation of A«had( — so) is reduced as 
So is decreased in comparison to that for s = ±M§ . This reduction may be anticipated, since 
from the form of (^) we see that the error mainly arises from uncertainties in the data for R(s') 
with s' < \sq\. 

Let us illustrate this point in more detail. If we compare the calculation of Aa had (— s ) data 
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y/s' interval (GeV) 


a = -6 GeV 2 


s = -15 GeV 2 


s = -50 GeV 2 


s = -Mf 


s = M| 


2m^ -0.81 


23.40 ± 0.48 


24.51 ±0.50 


25.07 ±0.51 


25.31 ±0.52 


25.31 ±0.52 


0.81 - 1.46 


11.31 ±0.50 


12.49 ±0.56 


13.14 ±0.59 


13.45 ±0.60 


13.45 ±0.61 


1.46 — 1.9 


< 


\ 5.90 ±0.40 
[ 7.05 ±0.72 


< 


f 7.27 ±0.50 
[ 8.70 ±0.89 


f 8.17 ±0.56 
\ 9.77 ± 1.00 


f 8.62 ±0.60 
\ 10.31 ±1.06 


< 


\ 8.62 ±0.60 
[ 10.32 ±1.06 


1.9-2.8 




f 6.95 ±0.44 




\ 9.70 ±0.61 


f 11.93 ±0.75 


f 13.24 ±0.83 




\ 13.26 ±0.83 


< 


[ 7.28 ±0.44 


< 


[ 6.67 ±0.16 


\ 12.42 ±0.75 


\ 13.77 ±0.83 


< 


[ 13.79 ±0.83 


2.8 - 3.74 


3.53 ±0.01 


5.68 ±0.02 


7.98 ±0.03 


9.65 ±0.03 


9.67 ±0.03 


3.74 - 5 


3.65 ±0.09 


6.67 ±0.16 


10.92 ±0.26 


15.06 ±0.36 


15.13 ±0.36 


5 — oo 


6.11 ±0.02 


13.36 ±0.09 


31.38 ±0.18 


169.99 ±0.52 


170.23 ±0.53 


u), (ft, if)'s, T's 


10.60 ±0.27 


13.37 ±0.38 


16.17 ±0.49 


18.73 ±0.58 


18.79 ±0.58 


A< a (-s) ■ 10 4 


1 


71.45 ± 1.13 
72.93 ± 1.41 


( 


93.05 ± 1.41 
94.91 ± 1.70 


f 124.76 ±1.64 
\ 126.85 ±1.92 


f 274.05 ±1.86 
\ 276.27 ±2.12 


( 


274.46 ± 1.86 
276.69 ±2.12 



Table 4: Explicit breakdown of the contributions to Aa^fs = — sq) in the spacelike region for 6 GeV 2 < s < Mf. A earn we 
show alternative results for the energy intervals 1.46 < Vs 7 < 2.8 GeV 2 and the final sum, where the upper (lower) braced entry 
corresponds to the use of inclusive (exclusive) data. The perturbative contributions here were evaluated with all u, d, s, c and 
b flavours in their active domains, and five light quarks contributing internal loops at C(a|) and O(ag). The scale is taken as 
fi = 20 GeV, the c pole mass as 1.4 GeV and the b pole mass as 4.7 GeV. For convenience we show in the last column the direct 
evaluation of Section 2, except that here, for consistency with the space-like evaluations, we use a fixed QCD scale fx = 20 GeV and 
five light quarks in the internal loops. The individual errors are combined as in Table 2. 
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Figure 5: Our determinations of a (M|) compared to two recent other determinations 
as two much earlier evaluations |l5j, 0] • 



] , as well 



with the direct evaluation of Aah a d(— Mf ) data , then essentially we make the replacement 



M 2 z R(s') 
s' + Ml 



R(s') 



sqR(s') 

s' + So 



(8) 



in the integrand of (|j), where for simplicity we consider s' <C M§ . Then we add to Aa n ad( 



-s 



\data 



-Ml) 



A«had(— so)]^ CD , as in (|5|). That is, if we compare the analytic 



the QCD term [Ac* na d( 
continuation determination, (|5|), of Aah a d(— Mf ) with the direct determination AcVhad( 
then effectively we make the replacement 



-M|) data , 



R{s 



I \ data 



s R{s') data + s'R(s')^ CD 
s + s' 



(9) 



for s' M|. Thus for s' s we keep all the data, while if s' ~ s we use pQCD to replace 
about half of the data, and for s' ^> so we discard almost all the data in favour of pQCD. 
Thus the lower that we can take so, the smaller the data contribution, and hence the smaller 
its contribution to the error on Aahad(i^i)- 

However before we can take advantage of the reduction of the uncertainty associated with 
the data, we must consider the error in the perturbative QCD continuation from s = —Sq to 
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Figure 6: x 2 fit as a function of the standard model Higgs mass, Mh, to the latest compilation of 
electroweak data, obtained using the 'inclusive' and 'exclusive' determinations of Ao!u^ d (M|) • 10 4 of 
274.3 (continuous curve) and 276.5 (dashed curve). The shaded zone to the left illustrates the energy 
interval where the Higgs has been excluded by direct searches at LEP2. 



12 




0.5 1 5 10 50 100 0.5 1 5 10 50 100 



Figure 7: Figure illustrating the efficacy of the Pade interpolation technique through threshold for 
fi = 20 GeV and a charm mass of m c = 1.4 GeV as a generic example. The C(a|) contribution to the 
Adler D-function is shown as high and low energy expansions, for (i) the pseudo-Abelian contribution 
containing no internal loops, (ii) the non-Abelain contributions containing triple gluon vertices, (hi) 
the contribution corresponding to the radiation of an internal light quark loop from a massive external 
quark loop, and (iv) contribution corresponding to the radiation of an internal massive quark loop, of 
the same mass scale as a massive external quark loop. The Pade threshold interpolation (continuous 
curve) becomes indistinguishable from the mass expansions away from threshold. 
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s = — M§. That is the error on 



s Q CD 



8(s ) = Aa had (-M z ) - A« had (-s 



(10) 



M l . , dU(s') 



—4ira I ds' 

-sq ds 



where IT, the hadronic contribution to the photon vacuum polarisation amplitude, satisfies 

Aa had (s) = -47raReII(s), R{s) = 127rlmn(s). (11) 

The error associated with the remaining analytic continuation round the semicircle from s = 
— Mf to s = Mf is much smaller and may be neglected, see (|6]). = 



To evaluate 5(s ) of ([TOf) we use the known expression for Il(s) to 0(a s ). For the 0(1) 



and O(as) contributions we use the full analytic formula fL9 |, which includes the dependence 
on the quark masses. The 0(a 2 s ) contribution is evaluated in terms of the high (m 2 /s) and 
low energy (s/m 2 ) expansions |20) using a (4/4) Pade interpolation^] for s ~ Am 2 [21J and, 
finally, the massless quark limit of the 0(a 3 s ) contribution is used. The expressions are valid 
for fixed coupling a s (/z 2 ). In Table 5 we show the individual contributions to S(s ) for a choice 
m c = 1.4 GeV of the pole mass of the charm quark and \i — 20 GeV of the QCD scale. 
Unfortunately there are appreciable uncertainties in the perturbative QCD determination of 
<5(so) arising from the sensitivity to the values taken for m c (and m^) and the QCD scale \i. 
In addition, in a recent paper Chetyrkin et al. [22] have evaluated the m A /s 2 term in the 



0(a\) contribution to R(s'). Of course knowing just the first two terms [23], |22j in the m 2 /s 
expansion is not sufficient to calculate the 0{a z s ) heavy quark effect, which comes mainly from 
the threshold region. However knowledge of these terms enables us to estimate the typical size of 
the 0{q z s ) mass contribution to be of the order of (0.2 — 0.5) x 10~ 4 . In total, these 'theoretical' 
uncertainties in the QCD contribution to (^) are comparable with the error presented in Table 1 
for the direct evaluation of Aahad(^f)- We conclude that although the error on Ac^at^—Sci), 
with So = 6 GeV 2 , is considerably improved in comparison to that for the direct determination 
of A«had(^|), nevertheless the uncertainties in the analytic continuation from — sq to Mf 
means that the accuracy to which Aat,ad(Mf ) is known has not been improved by the analytic 
continuation approach. 



4 Resolution of the "inclusive-exclusive" ambiguity 

We have seen that analytic continuation does not appear to allow us to reduce the uncertainty 
in the determination of Aa ha d(M|). However if we turn the analysis around we have the 
possibility to 

6 An example of the power of the Pade interpolation is shown in Fig. 7. To calculate the (D(a 2 s ) contribution 
to S(sq) we perform the appropriate integration of the Pade interpolation over the interval s = — sq to s = — Af§ . 
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1.62 


1.08 


0.78 


0.44 


0.18 


-0.10 


o(<4) 




0.26 


0.14 


0.10 


0.06 


0.03 


0.01 


5 QCD (s ) ■ 10 4 


u, d, s, c 


189.87 


168.53 


155.76 


137.76 


119.28 


32.31 



Table 5: The individual contributions to 5 QCD (s ) = [Aa had (-Mf ) - Aa had (s )] to 0(a|) from 
the u, d, s and c flavours. Note that the QCD contributions in the earlier Table |] also include the b 
quark. 



(i) distinguish between the inclusive and exclusive data for R(s) for -y/s < 2 GeV, 

(ii) constrain the value of the charm mass m c . 

To do this we study the difference between the 'direct' prediction for Aah ad (Mf) (shown in 
the last column of Table 4) and the values obtained via the analytic continuation method of 
eqs. OH) and (||). Let us denote the difference of the two determinations by d(s ), that is 

d(s ) = A« had (Ml) ^ - Aa had (M§) _ . (12) 



anal. cont. from so 



A self-consistent analysis requires that d(s ) ~ for all values of s . Of course the perturbative 



so 
GeV 2 


inclusive dat 
fi=10 ^ = 20 


a 

li = 50 


exclusive data 
[m = 10 /i = 20 n = 50 


6 


-0.31 


-0.12 


0.03 


0.43 


0.62 


0.77 


15 


-0.16 


-0.03 


0.07 


0.20 


0.33 


0.43 


25 


-0.13 


-0.01 


0.06 


0.11 


0.23 


0.30 


50 


-0.11 


-0.03 


0.04 


0.02 


0.10 


0.17 


100 


-0.07 


0.00 


0.04 


-0.01 


0.06 


0.10 



Table 6: The discrepancy d(so) x 10 4 of (|i"2| ) for space-like evaluations at s = —sq for three different 
scales [i (in GeV). In the first half of the table the inclusive data for R(s') is used in the region 
y/s* < 2 GeV, whereas in the second half the exclusive data are taken. 



QCD contribution depends on the value taken for the charm mass m c and the scale /i. We 

15 



therefore proceed in stages. First we remove the dependence on m c (and m&). We include 
only contributions from u, d and s quarks, and substitute for the data and resonances in the 
charm (and bottom) threshold regions with the values obtained from three-flavour perturbative 
QCD. The results for the discrepancy d(so) are shown in Table 6 for three different choices of 
the scale \i. It is immediately that, in general, if we use the inclusive R(s') data in the region 
x/s 1 < 2 GeV we obtain better agreement (that is a smaller discrepancy d(so)) than if we use the 
exclusive data. Moreover the scale fj 2 should be representative of the interval of continuation 
from s = —sq to s — — M§, and \x = 20 GeV is a reasonable choice. If we assume that the 
systematic discrepancy d(so) comes from a local region s' ~ s p then the additional contribution 
to the dispersion integral may be approximated by|] 

d{So) * £ / dS ' ' 6{S ' '- Sp) J^fsV) 

~ aR P (13) 
3tt(s + s p )' 

In fact the differences d(so) for the exclusive data at // = 20 GeV are well described by this 
simple pole form with 

y/s; = 2.1 GeV, R p = 0.8 GeV 2 . (14) 

This is consistent with the exclusive contribution being too large in the region y/s' ~ 2 GeV. 
We may conclude the three-flavour analysis of this section favours the inclusive data for R(s') 
for y/s' < 2 GeV and, moreover, gives a remarkably consistent description with d(so) ~ for 
different choices of Sq for scale choices in the region 20-50 GeV. 



5 Implications for the charm mass 

We now extend the 'discrepancy' analysis of the previous section to four-flavours and reinstate 
the data for R(s') in the charm threshold region (that is the J/ip, ip' and 3.74 < y/P < 5 GeV). 
We show the results in Table 7 for a range of choices of the charm mass m c , taking the scale 
H = 20 GeV. We see a systematic trend of the behaviour of d(s ) with m c and that the choice 
m c = 1.40 GeV gives good consistency for all So if the inclusive data are used in the region 
y/s' < 2 GeV. The numbers in brackets in Table 7 correspond to using the exclusive data up 
to y/s' < 2 GeV. There is no choice of m c that gives the same consistency as for the inclusive 
data. The optimum value appears to be m c = 1.34 GeV. 

The discrepancies d(so) were fitted to the pole form flilf ), and the parameters (the residue R p 
and pole position s p ) are given in Table 8. Again we see the inclusive data select m c = 1.40 GeV 
(for fj, = 20 GeV) and that as we depart from this value the additional pole contribution is 
such as to compensate for the poorer choice of m c . For the exclusive data we confirm that the 

7 We may use the unsubtracted form of the dispersion integral since in a difference calculation the subtraction 
constant will cancel. 
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m c { GeV) 


a[So = 6 GeV j ■ 10 


a(15) • 10 


a(25) • 10 


a(50) ■ 10 


a(100j ■ 10 


1.46 


0.57 (1.31) 


0.36 (0.72) 


0.23 (0.47) 


0.10 (0.23) 


0.04 (0.10) 


1.44 


0.39 (1.13) 


0.24 (0.60) 


0.16 (0.40) 


0.06 (0.19) 


0.03 (0.09) 


1.42 


0.20 (0.93) 


0.12 (0.48) 


0.07 (0.31) 


0.02 (0.15) 


0.00 (0.06) 


1.40 


0.00 (0.75) 


0.01 (0.37) 


0.00 (0.24) 


-0.03 (0.10) 


-0.03 (0.03) 


1.38 


-0.17 (0.57) 


-0.09 (0.27) 


-0.08 (0.16) 


-0.07 (0.06) 


-0.04 (0.02) 


1.36 


-0.37 (0.37) 


-0.21 (0.15) 


-0.16 (0.08) 


-0.11 (0.02) 


-0.06 (-0.01) 


1.34 


-0.57 (0.17) 


-0.33 (0.03) 


-0.24 (-0.00) 


-0.16 (-0.03) 


-0.09 (-0.03) 


1.32 


-0.74 (-0.01) 


-0.43 (-0.07) 


-0.31 (-0.07) 


-0.19 (-0.06) 


-0.11 (-0.05) 


1.30 


-0.94 (-0.20) 


-0.54 (-0.18) 


-0.39 (-0.15) 


-0.25 (-0.12) 


-0.14 (-0.08) 


1.28 


-1.11 (-0.38) 


-0.64 (-0.28) 


-0.45 (-0.21) 


-0.28 (-0.15) 


-0.15 (-0.09) 


1.26 


-1.32 (-0.58) 


-0.76 (-0.40) 


-0.53 (-0.29) 


-0.32 (-0.19) 


-0.18 (-0.12) 



Table 7: The discrepancy d(so) = 5 data (so) — (5^ CD (so) for a spectrum of charm pole masses and the 
lower QCD scale fi = 20 GeV. The entries (bracketed) correspond to the use of the interpolations of 
the inclusive (exclusive) R(s') data in the region \fs' < 2 GeV of the dispersion integral (||). 

value m c = 1.34 GeV is optimum, but that the pole compensation for other choices of m c is 
more more erratic. We repeated the whole analysis for scale \i = 50 GeV. The pole parameters 
which fit the discrepancy d(so) in this case are also shown in Table 8 (in the last two columns). 
For this choice of fi, the inclusive data give m c = 1.33 GeV whereas the exclusive data select 
m r = 1.26 GeV. 



Our determinations of m c refer to the pole mass of the charm quark. However the PDG [|12 
gives the value of the charm mass m c (fi = m c ) in the MS scheme, that is the running mass at 
scale m c . They quote m c (m c ) = 1.25 ±0.10 GeV, which is determined from charmonium and D 
meson masses. In our calculation the pole mass naturally occurs in the space-like continuation, 
with the 'running' included in the expression for the vacuum polarisation. The PDG value 
corresponds to a pole mass m c = 1.46 ± 0.11 GeV. We summarize the determinations^ in 
Table 9. 

Again we see that the results favour the inclusive measurement of R{s) in the region 
\fs < 2 GeV. First, the inclusive data satisfy the self-consistency check d(s ) ~ for different 
s for some value of m c , better than the exclusive data, see Table 7. Second, the prediction for 
the pole mass m c = 1.33 — 1.40 GeV is in better agreement with PDG expectations than our 
prediction m c = 1.26 — 1.34 GeV obtained using the exclusive data. 



8 Some years ago an analysis |24| of the moments of R c (s), obtained from e + e" 
m c = 1.34 ±0.02 GeV. 



annihilation, gave 
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'Inclusive' 


'Exclusive' 


'Inclusive' 


'Exclusive' 


m c { Lrev ) 


= zU LreV 




H = oil Lrev 


H = ou Lrev 


1.46 


(0.9,6.4) 


(1.8,4.4) 


(2.3,9.4) 


(3.1,6.6) 


1.44 


(0.6,5.8) 


(1.5,3.9) 


(2.0, 10.0) 


(2.8,6.6) 


1.42 


(0.3,4.5) 


(1.1,3.3) 


(1.7, 10.7) 


(2.5,6.5) 


1.40 


(0.0, -4.6) 


(0.8,2.4) 


(1.3, 10.8) 


(2.1,6.1) 


1.38 


(-0.4, 14.4) 


(0.5, 1.3) 


(1.0, 12.3) 


(1.8, 5.8) 


1.36 


(-0.7,9.8) 


(0.3, -0.8) 


(0.7, 17.8) 


(1.4,5.6) 


1.34 


(—1.1, 8.9) 


(—0.0, —5.1) 


(0.85, 117.3) 


(1.1, 5.3) 


1.32 


(-1.3,7.7) 


(**, **) 


(-0.1,-1.8) 


(0.8,4.6) 


1.30 


(-1.7,8.0) 


(-1.6,57.5) 


(-0.4,2.9) 


(0.5,2.9) 


1.28 


(-1.9,7.2) 


(-1.3,21.0) 


(-0.7,3.8) 


(0.2,0.8) 


1.26 


(-2.2,7.1) 


(-1.5,14.7) 


(-0.9,4.0) 


(-0.0,-8.3) 



Table 8: The parameters (R p , s p ) describing the simple pole fits, Eqn. (|i~3|) to the residual function 
d(so) for the spectrum of charm masses. The entry denoted by (**,**) corresponds to a residual 
sufficiently close to for all so to render the fitting procedure inappropriate. 

6 Summary 

Traditionally the value of the QED coupling on the Z pole has been determined by evaluating 
the dispersion relation, (f|), for Aah a d(s) at s — M§. In Section 2 we presented an updated 
calculation of Aahad(^l) using the latest available data for R = cr(e + e~ — > hadrons)/o"(e + e~ — > 
/i + yU~). The main uncertainty is the input for R(s') in the region 1.5 ^ vF ^ 3 GeV. The new 
(preliminary) BES-II data [0 have improved the knowledge of R(s') in the upper part of this 
region, so that the error on Aahad(Af§) is about ±2 x 10~ 4 corresponding to about ±0.03 on 
a _1 (M|) ||. However this error does not take full account of the effects of the discrepancy 



Source 


m c (GeV) 


inclusive 
exclusive 


1.33-1.40 
1.26-1.34 


PDG 


1.46±0.11 



Table 9: The pole mass of the charm quark determined from demanding self- consistency of the 
spacedike evaluation of A«h a d (that is requiring the discrepancy d(so) ~ for all so), compared to 
the PDG value | |12| |. Inclusive (exclusive) mean that R(s') is determined from inclusive data (sum of 
the exclusive channels) in the region yV 2 GeV. In both cases the lower and upper values quoted 
for m c correspond to scale choices fi = 50 and 20 GeV respectively. 
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between the inclusive measurement of R(s') and the sum of the exclusive channels in the energy 
region y/s' < 2 GeV. This discrepancy in R(s') leads, on its own, to a difference of 2.3 x 10~ 4 
in the value of Acthad(^i); see Table 2. Clearly it is important to resolve the dilemma. 

We confirm the general conclusion of Jegerlehner [fj that analytic continuation does not 
improve the accuracy of the determination Aahad(^f)- We find the evaluation of Aa na d(s) at 
the space-like value s = —s = —6 GeV 2 has a reduced error of ±1.4 x 10~ 4 . However the 
reduction in the error is more than offset by the uncertainty in the perturbative QCD analytic 
continuation from s = — sq to s = — M|, which arises from its dependence on the choice of 
charm mass and of the QCD scale. 

On the other hand the evaluation of @ at different space-like values s = — So proves to be 
very informative. For each evaluation Act!h a d(so) at a different, but sufficiently large, so, we 
can analytically continue to s — — M|, and then around the semicircle in the complex plane to 
s = M§, using perturbative QCD. We can compare these determinations of Aa n ad(^I) with 
the traditional method of directly evaluating at s — M§. In fact we found it convenient to 
study the difference 

d(s ) = Aa had (Mf) - Aa had (M|) (15) 

direct anal. cont. trom so 

as a function of sq. A self-consistent analysis requires d(so) ~ for all choices of sq. 

Indeed we found that the study of d(so) sheds light on the 'inclusive' versus 'exclusive' 
data dilemma, and provides evidence in favour of the former. But first we noted that the 
perturbative QCD analytic continuation was sensitive to the pole mass m c of the charm quark, 
as well as to the QCD scale \i. To eliminate the dependence on m c (and m&) we evaluated d(s ) 
using the data for R(s') in the region 2m n < y/p < 2.8 GeV and three-flavour perturbative 
QCD elsewhere. We performed the analysis using first the inclusive, and then the exclusive, 
data for y/s 1 < 2 GeV; in each case for three choices of the QCD scale. We found the 'inclusive' 
d(so) values were more self-consistent than the 'exclusive' behaviour of d(so). 

We exploited the sensitivity of the d(so) analysis to the pole mass of the charm quark in 
order to determine the value of m c . To do this we repeated the above procedure with the 
charm data reinstated and used four-flavour QCD. If the 'inclusive' data are used, we found 
that indeed there is a unique value of m c for which we obtain the same Aahad(^i) f° r the 
different space-like s = —sq values and for the direct evaluation at s = Mf . In this way, we 
determine the pole mass to be 

m c = 1.33 - 1.40 GeV, (16) 

if the QCD scale is \i = 50 or 20 GeV respectively. Just as in the three-flavour study, we found 
that the four-flavour analysis is less consistent if the 'exclusive' data choice is employed. 

In summary, we have presented quite a body of evidence to show that self-consistency of the 
results for the space-like and time-like evaluation of dispersion relation (ffl) selects the inclusive 
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measurements of R(s') in the region 1.46 < ^fs* < 1.9 GeV, as compared to the values of R(s') 
deduced from the sum of the exclusive channels. Thus we conclude that 

Aa^Mf) = (274.26 ±1.90) x KT 4 , (17) 

and hence that 

a'\M 2 z ) = 128.972 ± 0.026. (18) 

The corresponding results using the exclusive data, which are not favoured, are (276.49±2.14) x 
10 -4 and 128.941 ± 0.029. Precise measurements of R(s') in the energy region \fs' < 2 GeV 
are necessary to confirm our conclusion and, more important, to improve the precision in the 
determination of the QED coupling on the Z pole. 
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Note added in Proof 



The final BES measurements have just become available, see Ref. |[25|| . The measurements of 
R at a/s = 2,2.2,2.4 and 2.5 GeV are slightly higher than the preliminary measurements @. 
In fact the latter three points now lie on our input curve for R that is shown in Fig. 1. The 
point at v/i = 2 GeV has increased by about 5% to R = 2.18 ± 0.07 ± 0.18 [||. These small 
changes do not affect the results presented in this paper. 
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